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Magnetic relaxation measurements on a series of Yi_2:PriBa2Cu307_i {x — 0.13, 0.34, 0.47) sin- 
gle crystals were performed over a large field and temperature range in order to investigate the 
characteristics of the vortex matter across the second magnetization peak (SMP). The magnitude 
of the SMP varies non-monotonically with the Pr concentration; i.e., the irreversible magnetization 
normalized by its value at the onset field Hon displays a maximum for the x = 0.34 single crystal. 
The two characteristic fields, Hon and Hsp, follow different temperature T dependences: Hon oc T''°" 
and Hsp oc [1 — {T /Tc)^]'^"'' . The extracted values of the apparent activation energy U* and the 
creep exponent fi display a maximum at a field Hon < H* < Hsp. Their field dependences point 
toward the coexistence of both elastic and plastic creep for H > Hon- The degree of participation 
of each creep mechanism is determined by the charge carrier density, which controls both the elastic 
properties of the vortex matter and the pinning potential. 

PACS numbers: 72.20.My, 75.30.Vn 



A. Introduction 

One of the most challenging manifestations of the vor- 
tex matter in cuprate superconductors is the second mag- 
netization peak (SMP) or fishtail displayed in the isother- 
mal magnetization at an intermediate magnetic field H^p 
much lower than the irreversibility field ^i^i^iii^i^iL^i^ i ^ ° i ^ ^ 
Theoretically, it has been shown that at this field the 
flux line lattice crosses over from a dislocation-free Bragg 
glass structure present at low magnetic fields, which has 
algebraic decay of translational correlations, to a vortex 
glass ji^iii^ii^ The existence of this order-disorder transi- 
tion was later observed experimentally. '^'^^'A- The mag- 
netic field increases the entropy of the flux line system^^ 
beyond that introduced by the quenched disorder. Ad- 
ditionally, the field triggers the competition between the 
elastic Uei and pinning Upm energies at constant tem- 
perature (see, for example, Ref. 18). At low fields, the 
weak point disorder is important and the flux lines form 
an elastic quasi-ordered dislocation-free lattice. As the 
field increases, the elastic screening lengthfi^ hence the 
interlayer interaction, decreases, weakening the elastic re- 
sponse of the lattice and facilitating the invasion of topo- 
logical defects (dislocations) . The break of the long range 
order allows a better matching of the vortex lines to the 
pinning potential and, consequently, enhances the irre- 
versibility present in magnetization. This effect starts at 
the onset field Hon- Above this field, the vortex matter, 
assumed to be a lattice, displays a granular structure, 
namely, a network of dislocations that separate grains of 
vortex matter that still retain a higher elasticityp2£ The 
dislocation lines obey a characteristic dynamics including 
pinning and depinning processesJi^ 

Although such a behavior is universal, each class of 
cuprates shows a specific behavior in which the second 
magnetization peak is more or less pronounced, depend- 



ing on the relative values of the elastic and pinning ener- 
gies, which, in turn, along with the characteristic lengths 
(elastic screening length Lq and Larkin length) depend 
on the superfluid density Us (through the magnetic pen- 
etration length A because the effective mass is almost 
independent) anisotropy 7, and the disorder parame- 
ter 6dis (Refs. 19,23-25). Usually, in the absence of a 
special treatment e.g., particle irradiation, weak disorder 
due to oxygen vacancies is always present in cuprates. In 
YBa2Cu307_5, charge underdoping is obtained by oxy- 
gen removal from the Cu(l)-0(1) chains. This procedure 
induces also disorder consisting in different sequences of 
full Cu-O(l) and empty Cu(l) chains and an increased 
occupation of 0(5) sitesi^^ A second consequence of the 
underdoping is the fast increase in the anisotropy. Thus, 
the three parameters rig, 7, and Sdis are interconnected in 
a complex way. This makes it difficult to systematically 
study the effect of one of these material parameters on 
the evolution of the vortex matter from a quasi-ordered 
Bragg glass to a disordered amorphous or liquid state 
(see, for example, Ref. 27). 

It is possible to reduce the number of the above men- 
tioned parameters from three to two by applying a doping 
method that practically does not introduce additional ef- 
fective disorder. This is the case of Pr substitution for 
Y in YBa2Cu307_5. The density of free charge carriers 
is extremely sensitive to Pr even at optimal oxygenation 
due to carrier localization in the Fehrenbacher-Rice band 
as a result of Pr-0 hybridization ;^^'^^ while the disorder 
created by the Pr ions is less effective to pinning com- 
pared with oxygen. "^"^ It is remarkable that the integrity 
of the chains, hence the orthorhombicity, is preserved in 
this material, notwithstanding the change in the charge 
carrier density. So, the main source of the weak pinning, 
i.e. oxygen disorder, is maintained almost constant. The 
only effect of Pr on the pinning parameter results from 
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its dependence on the magnetic penetration length.!-!^ 
Therefore, the Pr substitution for Y in YBa2Cu307_5 
facihtates the investigation of the effect of the super- 
fluid density on the evolution and crossovers of different 
regimes of the vortex matter. 

In this paper, we explore the evolution of the vortex 
matter state in the temperature and field range where 
the second magnetization peak SMP is present by study- 
ing the magnetization and magnetic relaxation of a series 
of Yi_2:Pr2:Ba2Cu307_a single crystals. Our study has 
shown that the main ingredient that controls the evolu- 
tion of the vortex matter through the different regimes is 
the charge carrier density. The SMP is first enhanced and 
then suppressed as x increases. The reason for this be- 
havior is the softening of the elastic moduli, which makes 
the vortex lattice less stable to defect invasion. 



B. Experimental Details 

A series of Yi_a;Pra;Ba2Cu307_5 {x = 0.13,0.34, and 
0.47) single crystals was chosen to perform extensive 
magnetization measurements. The three single crystals 
for which data are presented here are platelets of size 
1.25 X 0.95 X 0.12 mm^, 0.75 x 0.4 x 0.1 mm^, and 0.83 x 
0.47 X 0.03 mm'^ with critical temperatures Tc of 82, 50, 
and 34 K, respectively. Details of the crystal growth are 
reported elsewherei^ 

Magnetic field H dependent magnetization M and re- 
laxation measurements were performed at different tem- 
peratures T in the reduced temperature T/Tc range be- 
tween 0.2 and 0.9 by using a Quantum Design supercon- 
ducting quantum interference device magnetometer with 
the external magnetic field applied parallel to the c-axis 
of the single crystal. We used the persistent current mode 
with a scan length of 40 mm, which guarantees an excel- 
lent magnetic field homogeneity. The single crystals were 
cooled in zero field to the desired temperature and the 
whole M{H) loop was recorded in increasing and decreas- 
ing fields with H steps chosen to get the finest details in 
M{H). After performing a hysteresis loop at a given 
temperature, the sample was warmed up to T >> Tc 
and zero-field-cooled to the next set temperature. 

For all the single crystals studied, the demagnetization 
factor D, calculated from the initial slope of the virgin 
hysteresis curves, was found to be higher than 0.9, hence, 
the magnetic induction B = H + 4tt{1 — D)M k, H. 
Therefore, we have used the magnetic field H throughout 
this paper instead of the magnetic induction B. 

Magnetic relaxation was measured in the normal dc 
mode by monitoring the time decay of the magnetic mo- 
ment. For these measurements, we zero-field-cooled the 
single crystal to the desired temperature, the magnetic 
field was then ramped to the target value, and the mag- 
netization was recorded as a function of time t every 175 
sec for about two hours. 



C. Results and Discussion 

Figure 1 shows the magnetization loops M{H) for 
three Pr concentrations measured at T/Tc « 0.3. These 
M{H) curves are almost symmetric around the M — 
axis, i.e., the reversible magnetization is nearly zero, and 
they display a similar field dependence. Specifically, the 
absolute value of the magnetization decreases with in- 
creasing field beyond the full penetration field, reaches 
a minimum at the onset field Hon, then it increases 
again reaching a second magnetization peak (SMP) at 
Hsp, and, finally, decreases to zero with further increas- 
ing the magnetic field. Previously, both Hon and Hsp 
have been related with the order-disorder transition from 
the dislocation-free Bragg glass to the disordered vor- 
tex matter with a glassy structure. For example, several 
theoretical studies take the peak field Hgp as the order- 
disorder Iinefi^ii2ii^i2^ while in experimental studies of 
YBa2Cu307_5, the order-disorder crossover has been 
taken either as the kink observed for Hon < H < Hgp 
(Refs. 15,32) or at Hon (Refs. 3,33). 

The goal in the present work is to study the effect of 
the charge carrier density Ug on the evolution of the vor- 
tex matter in the temperature and field range where the 
second magnetization peak SMP is present. A specific 
vortex state at a fixed T and H is determined by the n^- 
dependent interplay between the pinning landscape and 
the parameters that govern the stability of the vortex 
lattice, namely, disorder, critical current density Jc, and 
elastic moduli. The study of Pr-doped YBa2Cu307_5 
insures that Pr doping does not change the disorder but 
only decreases Ug. Hence, by studying this system, one 
studies the above mentioned interplay as a function of rig 
only. The decrease of Jc with decreasing ris is reflected in 
the continuous decrease of the irreversible magnetization, 
including its value at the SMP, with decreasing n^, i.e., 
increasing x (see Fig. 1). Hence, a plot of Mirr/Mirr,om 
the irreversible magnetization Mirr = {M~ — M+)/2 
(M~ and M'^ are the magnetizations measured in de- 
creasing and increasing magnetic field, respectively) nor- 
malized by the value of the magnetization at the onset 
field Mirr^om vs H/Hsp eliminates the effect of Jc{ns)- 
Therefore, such a plot reflects, at constant H , only the 
rig-dependent interplay between the pinning landscape 
and the elastic moduli, while, at constant doping x, it 
refrects the interplay between the pinning landscape and 
the iJ-induced disorder. 

A plot of Mirr/Mirr,on VS H / H sp mcasurcd at a re- 
duced temperature T/Tc = 0.3 and three different Pr 
dopings is given in Fig. 2. Note that the normalized 
irreversible magnetization in the SMP region displays a 
conspicuous nonmonotonic behavior as a function of the 
charge carrier density, i.e., it increases with decreasing 
charge carrier density down to ~ 0.165 [x — 0.34) 
followed by a decrease for even lower (higher Pr con- 
centrations). In fact, our magnetic investigations of the 
X = 0.53 single crystals have shown that the SMP is ab- 
sent at this Pr concentrationi^ 
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The non-monotonic behavior of Mirr/Mirr,on vs 
H/Hgp can be explained as follows. The elastic mod- 
uli C44 and Cge decrease with decreasing charge carrier 
density Us (increasing x)^ The decrease of the elas- 
tic moduli reduces the energy scale of the dislocations^^ 
[Ed = (C44C66)"^^^6/47r, where b is the Burger's vector], 
making easier their generation and, subsequently, a bet- 
ter matching of the flux-line system to the rig dependent 
pinning landscape. Therefore, at high charge carrier den- 
sities (e.g., the X ^ 0.13 single crystal), the flux-hnes 
system is rather stiff since the creation of dislocations re- 
quires a rather high energy. Hence, the density of dislo- 
cations is rather low and their relative contribution to the 
total irreversibility, i.e., pinning of the flux lines, though 
important, is not substantial. At low charge carrier den- 
sities (e.g., the X ~ 0.47 single crystal), the flux-lines 
system is soft since dislocations are easily created (C44 
and Cee are small, hence Ed is low). Nevertheless, de- 
spite the increased plasticity of the flux-lines system due 
to the increased number of dislocations, the values of 
Mirr/Mirr^on are Smaller than the ones for the low Pr 
doping (high Us) case due to the net decrease in the 
strength of the ris-dependent pinning landscape. How- 
ever, there is an optimum charge carrier density (e.g., 
the X ~ 0.34 single crystal) for which the softening of 
the flux lattice due to the decrease in the elastic moduli 
allows the optimum matching of the flux-lines system to 
the pinning landscape, maximazing the irreversible mag- 
netization Mirr/Mirr,on- 

The H dependence of Mirr/Mirr,on for a constant x 
displays a broad maximum around the field Hsp corre- 
sponding to the second magnetization peak. This sug- 
ests that the gradual process of field driven disordering 
gives rise to a continuous increase in the density of dis- 
location loops. Hence, elastic and plastic behaviors co- 
exist in different degrees over a large field range starting 
at Hon (below which only elastic creep is present, Ref. 
3) and ending far above H^p, with the elastic behav- 
ior dominant for H < H^p and plastic behavior domi- 
nant for H > Hsp. This result is supported by STM 
investigations which show the presence of dislocations 
even near Hon and by neutron scattering data, which 
have shown surviving Bragg peaks far above Hp. (Ref. 
35). Note that numerical simulations have shown that 
homogeneous flux line domains, i.e., systems displaying 
elastic creep, survive up to 0.8 x Hc2 (Ref. 36). 

The overall behavior discussed above is present over 
the whole temperature range 0.3 < T/Tc < 0.9. As the 
temperature increases, both characteristic fields Hon and 
Hsp shift to lower values, but they follow different T de- 
pendences; namely, Hon T''^"" (see Inset to Fig. 3) 
while Hsp oc [1 - (T/Tc)2]'''p (see Fig. 3), with = 
1.9, 2.1, and 2.2, and i^sp = 2.07, 2.53, and 2.23 for the 
X = 0.13, 0.34, and 0.47 single crystal, respectively. A 
power law dependence of Hon{T) was also reported for 
YBa2Cu307_5 single crystals but with lower values of 
the exponent namely between 1.09 and 1.3 (Refs. 
32,33). 



The Hsp (T) of the present study does not show the up- 
turn at high temperatures reported for optimally doped 
YBa2Cu307_5 (Refs. 6,8,15,37). Nevertheless, previ- 
ous measurements have shown that the Hsp (T) behavior 
is extremely sensitive to small deviations from optimal 
doping^i^i and that the suppression of the high temper- 
ature upturn for YBa2Cu307_5 occurs for S > 0.06. On 
the other hand, disorder introduced by electron irradia- 
tion was found to decrease, but not suppress, the high 
temperature upturn of Hsp{T) even for a fluence as high 
as 2 X 10^^ cm~^ (Ref. 38). Therefore, we conclude that 
the reduction of the charge carrier density by Pr substi- 
tution for Y gives rise to the suppression of the upturn 
in Hsp{T) at high T, hence to the monotonic decrease of 
Hsp with increasing T. Specifically, the decrease in the 
superfluid density Us with increasing doping x at high 
T increases the magnetic penetration length A, hence 
decreases Hsp {Hsp oc 1/A'', Ref. 2). Notice that the 
X dependence of the exponent Vsp shows a similar non- 
monotonic trend as the SMP itself. 

We also performed magnetization relaxation measure- 
ments on Yi_a;Pra;Ba2Cu307_5 at different fields and 
temperatures around the SMP in order to further study 
the evolution of the vortex matter with the charge car- 
rier density. Figure 4 is a plot of the time t dependence 
of the irreversible magnetization M^^r normalized by the 
first measured magnetization Mirr{tb) (tb 250 sec) for 
the X = 0.34 single crystal measured at a reduced temper- 
ature T/Tc = 0.4 and different reduced magnetic fields 
H/Hsp both below and above Hsp. We obtained Mirr{t) 
by subtracting the reversible magnetization, extracted 
from the M{H) loops (see Fig. 1), from the measured 
magnetization. The data of Fig. 4 are representative for 
all the single crystals measured. It is salient that the 
data do not follow a logarithmic t dependence. Hence, 
we analyzed these data in the framework of the collective 
creep theory in which^S 
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where Uq is the effective pinning potential, to is a macro- 
scopic quantity depending on the sample size and it 
should not be confused with the actual microscopic at- 
tempt time)^ and ^ is the collective creep exponent. The 
normalized relaxation rate S is obtained from Eq. (1) as 



SiJ,T,H) 
where 
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M„r(t) d\n{t) U*{J,T,H)' 



U*{J,T,H) = t/o + MsTln - 

to 



(2) 



(3) 



is the apparent activation energy, which is larger than Uq 
due to current relaxation. Note that, for convinience, we 
define 5* as a positive quantity. One can determine exper- 
imentally U* from Eq. (2) with the normalized relaxation 
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rate S obtained from the data of Fig. 4. In doing so, the 
relaxation rate is normalized to the initial magnetization 
Mirr{tb) rather than the time dependent magnetization. 
Since the variation in Mirr during the relaxation mea- 
surement is small, the error introduced is also small. ■^^ 
Before discussing the physics that U* would reveal, we 
discuss next the relationship between this apparent ac- 
tivation energy, which is accessible experimentally, and 
the actual activation energy. 

The actual activation energy U is a rather complex 
quantity involving not only a term due to the microscopic 
interaction between the flux lines and the pinning centers 
Uint, but also an extrinsic dependence on the distribution 
of the critical current density4S Miu et al.^ inferred that 



u = Um{j) HJc/J)- 



(4) 



Then, the relationship between U* and ?7, hence J7mt, 
can be obtained through the general dependence U — 
kBTln{t/to) (Ref. 42) along with Eqs. (2) and (4) (see 
the Apendix for more details) as: 
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Therefore, U* is always an overestimate of Umt since 
dUint/dJ is negative. Nevertheless, Eq. (5) shows that 
the approximation of Uint with U* is valid as long as 
the current density is close to Jc and it breaks down for 
J <^ Jc i.e., for H > Hsp. Thus, a plot of U*{J) deter- 
mined from Eq. (2) with S obtained from the relaxation 
data in the regime where J is not too far from Jc, i.e., 
for not extremely long relaxation times and for magnetic 
fields around Hsp, gives an accurate information on Uint, 
hence, on the evolution of the vortex matter when the 
temperature and magnetic field are swept. 

The plot of U* vs Mirr {Mirr oc J) IS shown in Fig. 5 
for the X = 0.13 single crystal measured at T/Tc « 0.4 
and different values of the reduced field H/Hsp. It is 
salient the different evolution of the activation energy 
below and above Hsp. Specifically, for fields smaller than 
Hsp (open symbols), U*{J,H) increases rapidly as the 
current J (or equivalently Mirr) decreases, an expected 
behavior for an elastic vortex system in the collective pin- 
ning regime. Above Hsp (filled symbols), the increase of 
U* with decreasing current density becomes slower and 
slower suggesting a smooth crossover at Hsp to another 
regime, most likely a regime dominated by the fast pro- 
liferation of dislocations. 

We have observed a similar behavior of U*{J) for all Pr 
concentrations studied. However, for the same reduced 
field range and for T/Tc « 0.4, the range of values of 
U* decreases with increasing Pr concentration; e.g., U* 
varies between 200 - 1000 K for the x = 0.13, 100 - 600 
K for the x = 0.34, and 50 - 400 K for the x = 0.47 
single crystal. The decrease in the range of U* with 
increasing x is in agreement with previous results on 
Yi_2;Pr2;Ba2Cu307_5 polycrystalhne pellets.^ 

As a function of field, U* (or equivalently S, see Eq. 
(2)) systematically shows, at any relaxation time t, a 



maximum (minimum) at a field H*, which is between 
Hon and Hsp [see, for example, the Inset to Fig. 5 for 
S{H) for the x = 0.34 single crystal]. Therefore, the 
slowest magnetic relaxation takes place at a field value 
just below Hsp and not at Hsp. At fields higher than H*, 
S{H) increases almost linearly. 

The apparent activation energy U* is a nonlinear func- 
tion of time, but a good estimate of its field dependence 
in a certain time window ti to t2 can be obtained by re- 
placing the derivative in Eq. (2) with finite differences. 
This average apparent activation energy U is given by 



Tj*ru\ 1 T^l (^"^ \^irr{h, H) 

U iH) = knTH-)^^^-^ 



(6) 



where AMirr{H) ~ Mirr{H,ti) — Mirr{H,t2). Figure 
6 is a log- log plot of U (H), as extracted from relax- 
ation measurements over the time window 600 sec < t < 
4000 sec, for three Pr concentrations {x — 0.13, 0.34, and 
0.47). U (H) displays again a maximum at H* . The rel- 
ative difference {Hsp — H*)/Hsp between H* and Hsp 
is maximum for the x — 0.34 single crystal, for which 
the value of the second magnetization peak is enhanced 
(see Fig. 2). For H < H* , U increases with increas- 
ing magnetic field, which is consistent with the elastic 
(collective) creep mechanism. For H > H*, U {H) de- 
creases roughly as a power law with increasing H, i.e., 
U*{H) oc H-", with p = 0.6 and 0.4 for the x = 0.34 
and 0.47 single crystal, respectively, which indicates plas- 
tic vortex creep^ and a more abrupt drop for the x = 0.13 
single crystal. A simple model for plastic pinning^ yields 
an exponent v = 0.5, while other reports give v = 0.55 
for BizSraCaCusOg (Ref. 44), j/ = 0.7 for YBazCuaOT-i 
(Ref. 2) and HgBa2Cu04+5 (Ref. 45), and i/ « 0.9 for 
Tl2Ba2CaCu208 (Ref. 46). Hence, as expected, these 
U {H) data are consistent with the U*{J,H) data of 
Fig. 5, but, additionally, they give the quantitative de- 
pendence of the activation energy on the magnetic field. 
Based on these data, we also conclude that, at least in 
the case of Pr-doped YBa2Cu307_5, the order-disorder 
crossover is best given by H*{T) {Hon < H < Hsp), 
eventhough, as discussed above, the two characteristic 
fields Hon and Hsp are the ones which previously have 
been related with the order-disorder transition. 

Based on our above conclusion, leveling the elastic and 
plastic energies engaged in the equilibrium of the flux 
lines, i.e., Uei = Upi, at H* gives H* oc l/T^A-* oc 
[{1 - (r/Tc)4}/T]x, with X = 2. The fit of H* vs 
[1 - {T/Tc)'^]/T for 0.3 < T/Tc < 0.8 is shown as a log - 
log plot in the inset to Fig. 6. The values of the expo- 
nent X are 1.44, 1.27, and 1.54 for the x = 0.13, 0.34, and 
0.47 single crystal, respectively. The good fit of the data 
with the above expression supports our conclusion that 
the elastic to plastic crossover takes place at H*. It is 
interesting to note that the exponents are non-monotonic 
with increasing x, with the lowest value for the x = 0.34 
sample. 

Since, as shown above, the activation energy is a func- 
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tion of current density, temperature, and magnetic field 
[see Eq. (3)] , changes of any one of these parameters drive 
continuously the vortex matter into different elastic and 
plastic creep regimes. Hence, magnetic relaxation data 
give information about a specific flux-creep regime for 
a given T and H through the critical exponent /i(r, H) 
present in Eq. (1). A fit of the Mirr vs \x\t data with 
Eq. (1) for different T and H gives ii{T,H). 

Figure 7 and its insets are plots of /i(T, H) for differ- 
ent charge carrier densities, i.e., Pr doping. Note that 
fJ.{H) displays a peak at the same magnetic field value 
H* at which U (H) is maximum. Therefore, as dis- 
cussed above, that elastic pinning mechanism dominates 
for Hon < H < H* while the plastic mechanism domi- 
nates for H > H* . Also, note that /j, decreases with in- 
creasing temperature and decreasing charge carrier den- 
sity (increasing x). 

The collective (elastic) creep theory^'^ predicts that 
fi = 1/7 for single vortex creep (at high current and 
low field), = 5/2 for small vortex-bundles creep (at 
intermediate current and field), = 1 for the creep of 
intermediate vortex bundles, and = 7/9 for the creep 
of large vortex bundles (at low current and high field). 
As the data show, /i(^f) does not follow exactly these 
theoretical predictions in the regime where flux lines are 
expected to behave elastically, i.e., for Hon < H < H*. 
For example, n > 2 but smaller than 2.5 for the lowest 
measured temperature [T/Tc = 0.3) even in the case of 
the a; = 0.13 single crystal, which is expected to have the 
strongest elastic response. Hence, although these values 
indicate that the relaxation of the flux vortices is mainly 
due to the creep of small vortex bundles, the admixture 
of the plastic contribution limits to values smaller than 
2.5. Similar deviations from the theoretical exponents for 
plastic creep^ are clear for H > H*. 

The direct (inverse) correlation between the values of 
the creep exponent /it and the the charge carrier density 
(Pr doping) as well as the the decrease of the value of 
fi with increasing T (see Fig. 7 and its insets) reflects 
the decreased role of the vortex lattice properties over 
single vortex behavior with decreasing Us or increasing 
T as a result of the weakening of the elastic moduli with 
decreasing Ug [Cge oc oc ns{x)] or increasing T. As a 
consequence, the crossover magnetic field Hsb from sin- 
gle vortex to small bundle collective pinning increases 
with decreasing (increasing x) or increasing T since 

Hsb oc A®/-^ oc nj^^^ (Ref. 23). Hence, H^b cannot be 
reached in strongly under-doped single crystals before 
the full crossover to plastic pinning. This explains the 
decrease of the value of the exponent fi at H* in Fig. 7 
and its insets with increasing x oi T. 

For magnetic fields higher than H* , decreases mono- 
tonically with increasing H. However, the value of 
/i = 10/21, representative for the plastic creep of the 
lattice, is accessible only at high temperatures. At lower 
T one needs magnetic fields higher than the one available 
(5 T) in order to be able to detect this regime. Addition- 
ally, the elastic contributions do not vanish completely 



in the plastic regimsi^ 



D. Summary 



In summary, we investigated the evolution of the sec- 
ond magnetization peak (SMP) with the charge carrier 
density by an appropriate doping that avoids the change 
of the quenched disorder in a sensitive way. For this goal, 
we carried out magnetization and magnetic relaxation 
measurements on a series of Yi_a;Pra;Ba2Cu307_5 single 
crystals in which the concentration of Pr ions controls 
the charge carrier density ns{x). We have found that 
the quenched disorder is necessary for the existence 
of irreversibility and of the SMP, but the principal 
ingredient which controls the evolution of the vortex 
matter through different regimes is the charge carrier 
density. Specifically, we have found that the SMP is 
broad and its magnitude is non monotonic with the 
amount of doping: it increases with decreasing the 
charge carrier density up to a doping around x = 0.34 
followed by a decrease with further decreasing 
(increasing Pr concentration). The two characteristic 
magnetic fields, the onset field Hon and the field Hgp 
corresponding to the SMP decrease with increasing 
temperature T, but they follow different T dependences: 
Hon oc T""^ while Hsp oc [1 - {T /T^)^]""- , with the 
exponent v^p following the same nonmonotonic trend 
as a function of tIj as H^p. Within the collective creep 
theory, we determined the apparent activation energy. 
Its evolution with J has shown that the vortex system is 
predominantly elastically pinned below Hgp, while above 
Hgp there is a smooth crossover to a vortex regime most 
likely dominated by the proliferation of dislocations. 
The field dependence of the average apparent pinning 
potential U displays a maximum at a magnetic field 
H* , with Hon < H* < Hsp, which is consistent with the 
presence of an elastic (collective) creep mechanism at 
low fields and plastic vortex creep at high H values. The 
transition from the Bragg glass to the dislocation rich 
vortex system occurs gradually and extends on a rather 
large field range. For this reason, we propose that the 
order-disorder line must be defined by the maximum of 
the average activation energy, which is located at H*, 
below Hsp but above the inflection point of the M{H) 
curves. 
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E. Appendix 

As mentioned in the main text, Miu et al^ inferred the 
following current dependence of the activation energy: 



U{J)^U,^tiJ) ln( 



Jc 



Then, 



dJ ^ dJ 



d(ln - In J) 



dJ 



J 



\dJ ) dJ \j) 



(7) 



As shown by Eq. (2), the normalized relaxation rate is 
1 dM„r{t) kgT 



Mirrit) dln(t) U*{J,T,H)' 



(8) 



Also, the actual activation energy given by 
U = kBTln{t/to) 

implies 

IhTt ^ 

Since Mi^r cx J, Eq. (8) becomes 

I dJ dU 1 



J dint dlntU*' 



hence, 



U* = ~J 



,dU 
dJ' 



Using Eq. (7), the above equation gives 
dU 



U* 



-J 



dl 



(9) 



(10) 
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I. FIGURE CAPTIONS 

Fig. 1 (Color online) Magnetic hysteresis loops for 

Yi_^Pr^Ba2Cu307-5 {x = 0.13,0.34 and 0.47) single 
crystals, measured at the same reduced temperature 
T/Tc w 0.3. The arrows indicate the position of the 
onset Hon and second magnetization peak Hsp fields. 



Fig. 2 (Color online) Plot of the irreversible mag- 
netization Mirr normalized to its value Mirr,on at the 
onset of the second magnetization peak, as a function 
of the reduced field H/Hgp for Yi_j:Prj:Ba2Cu307_5 
{x = 0.13,0.34 and 0.47) single crystals measured at the 
same reduced temperature T/Tj. k. 0.3. Inset: x depen- 
dence of the absolute value of the irreversible magnetiza- 
tion at the second magnetization peak Mirr,spi measured 
at the same reduced temperature. 

Fig. 3 (Color online) Log-log plot of the second mag- 
netization peak Hsp vs reduced temperature T/Tc for 
Yi_^Pr^Ba2Cu307_5 {x = 0.13,0.34 and 0.47) single 
crystals. Inset: Log-log plot of the onset field Hon vs 
T/Tc for the same single crystals. 

Fig. 4 (Color online) Logarithm of time t evo- 
lution of the irreversible magnetization Mi^r of an 
Yo.66Pi"o.34Ba2Cu307_5 single crystal measured at the 
reduced temperature T/Tc = 0.4 and different field values 
around the second magnetization peak. Open symbols 
below Hsp, and partially filled or close symbols above 
Hsp represents 1.06, 1.21, 1.51, 1.82, 2.12, 2.42, 2.73, 
2.03, 3.64 as shown by arrow. 

Fig. 5 (Color online) Plot of the apparent ac- 
tivation energy U* , obtained from relaxation mea- 
surements, vs irreversible magnetization Mirr of an 
Yo.87Pro.i3Ba2Cu307_5 single crystal measured at dif- 
ferent applied magnetic fields and at the reduced tem- 
perature T/Tc ~ 0.4. Empty and filled symbols are for 
H < Hsp andiJ > Hsp, respectively. Solid lines are guide 
for the eye. Inset: Field dependence of the relaxation rate 
S measured at a time t w 6000sec. 

Fig. 6 (Color online) Log-log plot of the field H 
dependence of the average activation energy U* of 
Yi_^Pr3;Ba2Cu307_i (x = 0.13,0.34 and 0.47) single 
crystals measured at a reduced temperature T/T^ = 0.4. 
The lines are fits of the data with a power law. Inset: 
Temperature dependence of H* (T) determined from the 
minima of S{H). 

Fig. 7 (Color online) Field H dependence of the 
relaxation exponent /x of an Yo.87Pro.i3Ba2Cu307_5 
single crystal mcasiircd at different reduced tempera- 
tures (Tc = 82 K). Insets: The same plot for the (a) 
Yo.66Pi'o.34Ba2Cu307_5 single crystal {Tc = 50 K) and 
(b) Yo.53Pro.47Ba2Cu307-5 single crystal (T^ = 34 K). 



